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Abstract 

Using Finite-Size Scaling techniques, we numerically show that the first irrele- 
vant operator of the lattice theory in three dimensions is (within errors) com- 
pletely decoupled at A = 1.0. This interesting result also holds in the Thermody- 
namical Limit, where the renormalized coupling constant shows an extraordinary 
reduction of the scaling-corrections when compared with the Ising model. It is 
argued that Finite-Size Scaling analysis can be a competitive method for finding 
improved actions. 
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1 Introduction 

The study of the scale-invariance of a system with many-degrees of freedom is of relevance 
both for the study of phase transitions in Condensed-Matter Physics and for High Energy 
Physics. The Renormalization Group (RG) is the central concept for these investigations. 
In both fields it can be turned onto a powerful tool when combined with Monte Carlo 
(MC) simulations In this kind of investigations one is only interested in universal 
properties (in RG language, the leading singularities described by the relevant operators). 
However, the action used in the simulation is also coupled to irrelevant operators. Their 
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effects can only be neglected in the limit of infinite lattice size and infinite correlation- 
length. This systematic error is at present the main difficulty to extract meaningful 
results from the simulation. For phase-transitions studies, the quantitative consideration 
of the scaling-corrections induced by the irrelevant operators in Finite-Size Scaling (FSS) 
studies was started in Ref. 0. At present, the problem is quite well understood, both 
below 0, 1^, and at the upper critical dimension 0. 

An action completely decoupled of the irrelevant operators is called a perfect action |§| . 
From the RG point of view, a perfect action is on the RG-trajectory that leaves the 
fixed point along the (supposedly unique) relevant direction (the so-called Renormalized 
Trajectory). It can be easily understood that the effort of finding such an action may be 
rewarding. In fact, the quest for perfect actions has been very active in the last years 
(see Ref. and references therein). The perfect action program may be summarized 
as follows. One first choose a real-space RG transformation, usually keeping a tunable 
parameter. A finite-dimensional coupling space is chosen, and the RG trajectory is tried 
to fit in it. The tunable parameter is played with, in order to minimize the truncation 
errors due to the finite-dimension of the coupling-space [ITU . 



The main difficulty in the perfect action program described above is related with the 
arbitrariness introduced by the RG transformation and by the coupling-space chosen to 
parameterize it. Indeed, given a reasonable coupling-space, it is expected to contain a 
set of points for which scaling-corrections are minimal (this is a working definition of a 
"perfect" action). It is clear that the most efficient way of locating it cannot be trying 
different RG group transformations until one of them falls just on it. We will argue that 
FSS analysis provides a simple and efficient way to locate this privileged set of couplings. 
To put the argument at work, we have investigated the lattice scalar theory in three 
dimensions. We will show that in this simple two dimensional coupling space it is possible 
to put the action's coupling to the first irrelevant operator below the statistical errors. 
However, the usual lattice gauge-theories simulations are not done in the FSS regime 
(correlation length ^ much larger than the lattice size L), although some interesting 
calculations have been performed For this reason, we have also investigated the 

scaling properties of our "optimal" action in the opposite regime, L S> ^. We have found 
again that the effects of the leading irrelevant operator lies below the statistical errors. 
As a byproduct, we are able to reconcile some discrepancies between MC calculations of 
the renormalized coupling constant at zero external momentum |13[, with a previous 
field-theoretic estimate |jl4 . 



Very recently a similar MC study has been published |15| for the = 1 Ising model. 
The freedom introduced is quite close to ours (we let the spin to be a real number, they 
only allow it to be 1,0 or —1). Thus, it is worthwhile to briefly comment on the differences 
between the two works. To start with, they use as input the value of a numerically 
determined universal quantity (basically a Binder cumulant), in order to numerically 
tune the action. On the contrary, we do not require a previous knowledge of the Binder 
cumulant. Next, they use their improved action to produce a strong error reduction on the 
estimates of the Ising critical exponents. We believe this error reduction to be unjustified 
(see Ref. 0, and section]^). 
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2 Finite Size Scaling and the best action 



To frame the discussion, let us consider a model with only two relevant parameters, the 
"thermal field", t, and the "magnetic" field, h, like for instance the Ising model. The 
free-energy of a sample of linear size L can be written as P, ffl 



fit, h, {u,}, L-') = g{t, h, {u,}) + s-<^/,i,g(s^*t, sy^h, {ujsy^}, s/L), j > 3. (1) 

where s is the scale of an (unspecified) RG transformation, and fsmg stands for the sin- 
gular part, g being analytical. The y's are the eigenvalues of the RG transformation, so 
that yt,yh > 0, but > 2/3 > ?/4 > . . . . A basic assumption of this approach is that the 
RG-evolution of the couplings is the same as in an infinite system. It is quite common to 
write z/ = 1/yi, r] = 2 + d — 2yh and 00 = —y^. From the free-energy, the thermodynam- 
ically interesting observables follow by taking derivatives. Let C,{L,t) be any reasonable 
definition of the correlation length in a finite lattice (see for instance Ref. |T6|). One can 
obtain a general expression for any observable, O diverging in the thermodynamical limit 
like t-^o; 



0{L,t) = L^o/^ 



L 



(2) 



where ^00 is the correlation-length in an infinite lattice. Eq. (^ is most interesting deep 
in the scaling region (^00 ^ L), where the 0{^^) corrections can be safely neglected. 
It is clear, that in Eq. (^ we have only kept the corrections due to the first irrelevant 
field, M3, (in fact, a full series in L~'^ is to be expected) but similar corrections arise from 
the others. A radically different correction is produced by the analytical part of the free 
energy, g. When one takes derivatives in Eq. (|ID with respect to the "magnetic" field, a 
not diverging term follows from g. This implies that for several important operators like 
the susceptibility, the Binder cumulant, or the correlation length, a correction like L~"'^'^ 
should be added in Eq. (H), if the multiplicative structure is to be kept. 

One can exploit Eq. (0), by comparing the measures taken in two lattices 



Q 



o 



{0{sL,t)) 
{0{L,t))- 



(3) 



at the value of the "temperature" , for which the correlation length in units of the lattice 
size is the same in both 



o 



(4) 



where Aq^ is a constant. In this way, not only the critical exponent xq/i^ can be measured 
but also information on the scaling corrections Aq^L~'^ can be extracted 0, ^ |l^ 18, T9[ 



Here we encounter an objective property that a perfect action should have: Aq^ should 
be zero for all the observables. Notice that this does not mean that there are no scaling 
corrections. There is no way to avoid the 0{L~^^'^) 0{L~'^) induced by the analytical 
part of the free energy (they will be there, even if all the irrelevant fields in Eq. (|l|), uj 
are set to zero). But we are still not done, because to efficiently locate the points for 
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which Aq^ = something should be known about xq/i'- The most simple procedure is 
to choose for O a quantity having xo = 0, like a Binder cumulant, B (see section ^ for its 
definition). 

Therefore, having a simple model, like for instance the Ising model, a simple strategy 
for improving its scaling behavior is the following. The original model has only a tunable 
coupling /3, that reaches a critical coupling, [3c (in Eq. (|l|) t can be identified with the 
"reduced temperature" {P — Pc) / Pc) ■ We need to extend the action with another coupling, 
A, which can be for instance a next-to-nearest neighbors coupling, a X(j)^ term, etc. The 
critical point Pc will extend into a critical line in the (/?,A) plane, /3c(A). If there is 
something like a RG fixed point in the (/3,A) plane, one should have 

^ B{P,X,sL) ^ aP,X,sL) 
B{P,X,L)'''' e(AA,L)' 

which is Eq. (^) in the absence of scaling corrections. We shall show that we can tune A 
to the condition given in Eq. (^), with moderate numerical effort, and that this yields a 
quite strong reduction of scaling corrections. 

One could object that there is no reason for the Fixed-Point to lie in such a small 
coupling-space. The answer is twofold. First, there is quite clear empirical evidence 
that this indeed happensQ for instance in the two dimensional Ising model |2^, the 0(4) 
Non Linear a Model in three dimensions and in the three dimensional site diluted 
Ising model, at 80% spin-concentration [0. Second and most important, we are only 
demanding that there exists some RG transformation whose Fixed-Point is close (in some 
sense) to the (/5,A) plane, but we are not burdening ourselves with the task of finding it. 
If it exists, we should be able of finding the fixed-point coordinates with Eq. (^. And 
that can be done with rather small lattices. 

3 The model 

We have considered the scalar A0^ theory in a three dimensional cubic lattice of side L, 
with periodic boundary conditions. The action is given by 

s = pj2 + + - 1)' ' (6) 

<i,j> i i 

where the 0's are real variables. The A — *• cxo limit of the model (|^) is the Ising model 
for magnetism. At fixed A, for low values of P the system is in a Z2-symmetric (param- 
agnetic) phase. There is a critical line, PdX), separating the paramagnetic phase from 
the ferromagnetic phase. The full critical line is second order, and it belongs to the Ising 
Universality Class, excepting the A = end-point which is the Gaussian model. The MC 
simulation of (^ is greatly eased by the cluster-methods PT|] . 

^From the numerical point of view, a truly perfect action and an action decoupled from the first 
irrelevant operator cannot be really distinguished, since one anyway has tiny but measurable analytical 
corrections. 
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The observables to be measured are easily defined in terms of tlie Fourier transform 
of the spin-field 

0(fc) = ^E<^(^)^''^ (7) 

X 

where V = Li^ is the lattice volume. We have the susceptibilities at zero (x) and minimal 
momentum (F): 

x = v {\m\') ^F = \ Yl (iSfc)i' 

iifcii=¥ 



and the finite-lattice correlation-length |T6| and Binder cumulant. 

The renormalized coupling constant at zero external momentum can be readily calcu- 

(10) 



lated lit 



(7R, = 2S(L/0^ , D = 3. 
We also measure the nearest-neighbor energy 



E 



"J / ' 



in order to calculate /3-derivatives and extrapolations |^. Let us summarize the critical 
behavior of the different observables: 



7 = v{2-7]) 



X 



, Xg^^ = 1 + 

the critical exponents being 

u = 0.6294(10) , 7] = 0.0374(12) , to = 0.87(9). 
These MC results are in excellent agreement with the last series estimates 



(12) 



(13) 



A final remark is that 



9r = lim lim gR{L,(3) ^ lim lim 9^(1,(3) = g^. 



(14) 



as Eq. (H) shows. 
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4 Optimizing A 



For a first estimate of tlie Fixed-Point location, we liave simulated lattices L = 4, 8 and 
16 in a mesh of A values (A = 0.25, 0.5, 1.0, 1.5 and 2). Results in the Ising limit were also 
available from previous work 0. We have measured Qb\q^=2 lattice-pairs (4,8) and 
(8,16), and check for the condition in Eq. (^ . Given the strong statistical correlation 
between and Qb, this quotient can be very precisely calculated. One can obtain a 
0.2% accuracy in the (8,16) pair, with nine Pentium Pro hours. The original idea was to 
use the mesh as starting point for a bisection search of the optimal A, but the A = 1 data 
fulfilled the condition in Eq. (^) to the achieved accuracy. Moreover, for the lattice pair 
(8,16), at A = 1.5, Qb was at three standard deviations below 1, and at A = 0.5 was five 
standard deviations above. From this, one can estimate the optimal A to be 1.0(1). The 
total CPU time was less than two hours of the 32 Pentium Pro machine RTNN. 

A more precise location of the optimal A would require not only much more statistics, 
but also to simulate larger lattices to avoid higher-order scaling corrections (if needed, 
one could resort to a reweighting method in A). Nevertheless our main scope, has been 
to check if a zero value of Aq^ (see Eq. (^)) does imply the vanishing of the 0{L~^) 
corrections for the other quantities of interest. In other words, if the first irrelevant 
operator is decoupled. For this, we have simulated at A = 1, and also at A = 0.5 and 
A = 0.25 for comparison. The total CPU time employed has been about three RTNN 
months. We remark that such a time-consuming investigation is not an essential part of 
an action improvement program, but rather a consistency check. 

In practice, one is mainly interested in the thermodynamical limit, where the mea- 
sures are quite less accurate than in the FSS region. For instance one could be interested 
in measures of exponential tails of propagators, or in measures of renormalized coupling 
constants. In the model considered in this work maybe the more interesting universal 
quantity is g^. Even with the high accuracy allowed by the improved cluster estima- 
tors we will show that a 10% error in the optimal A is enough to cancel the leading 
scaling-corrections. In fact the ^ = 2.5 data are in the continuum limit to our accuracy 
(2%). In this regime, we have used 6 RTNN days. 

5 Numerical Results when ^ L 

We have simulated in lattices L = 8, 12, 16, 24, 32 and 64 for A = 1.0, 0.5 and 0.25. The 
Elementary Monte Carlo Step (EMCS) was composed of 10 single-cluster flips, followed 
by a full-lattice Metropolis sweep. We have taken 10'' measures, separated by four EMCS, 
excepting L = 64 where only 2.5 x 10® measures were taken. Our values for Qb\q^=2 
quoted in table |l|. From Eq. (^) we can obtain an (L dependent) estimate of u and rj, and 
also of B{L)\q^^2 |;|q^_2- The results are displayed in Fig. ^ where a quadratic fit 
in L~'^ is presented. The results for the linear term in L^'^ of the A = 1 data are (fitting 
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Figure 1: Universal quantities as a function of L~'^ [u — 0.87), for the four models 
considered. The lines correspond to quadratic fits, constrained to yield the same infinite- 
volume extrapolation. 
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A = 0.25 


8 
12 

16 
24 
32 
48 

61 


0.98146(20) 
0.98703(23) 

0.99038(23) 
0.99324(21) 
0.99476(28) 
0.99612(28) 

0.99702(31) 


1.00094(25) 
1.00075(24) 
1.00062(26) 

0.9997(5) 


1.01372(30) 
1.01088(32) 
1.00817(29) 

1.0045(5) 


1.02949(40) 
1.02251(35) 
1.01840(38) 

1.0095(7) 



Table 1: Quotient of the Binder cumulant measured in lattice pairs {L,2L) where the 
correlation-lengths in units of the lattice size are the same. Scaling corrections appear as 
a value different from one. 
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Figure 2: Universal quantities as a function of L^"'^", for the A = 1 data. The hnear 
behavior shows that there is no evidence for scahng- correct ions others than the analytical 
ones. 



independently from the other A values) 

Aq^^^ = 0.004(44), Aq^ = 0.016(25), 

%L = 0.022(23) , Ab = 0.010(21) . ^ ' 

Thus, it seems quite clear that for the four quantities shown in Fig. |^, the coefficient of the 
0{L~^) corrections changes sign close to A = 1.0 and that the A = 1.0 data are compatible 
with a "perfect" behavior. That would mean that the scaling-corrections may arise from 
other irrelevant operator, and from the analytical corrections 0{L~'^^'^). Plotting our data 
against the analytical corrections, a nice linear behavior is found (see Fig. |^). Therefore, to 
our precision A = 1 seems perfectly decoupled from the first irrelevant operator. Operators 
with uo' > uo cannot be resolved from the analytical corrections. At this point, one could be 
tempted of using the fits of Fig. ^ to reduce the error in the infinite volume extrapolation 
for the critical exponents, as it is done for instance in Ref. [|r^. Followed blindly, this 
criterion reduces the errors in a factor four, or larger when compared to the fit in Fig. [l|. 
But we believe this reduction to be unjustified, because we have tuned A to Eq. (||) 
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numerically, so the condition Aq^ = can only be expected to hold within errors, as 
shown in Eq. (|15[). Therefore, the more conservative estimate given in Eq. (|13D is to be 
preferred. If a procedure were designed (numerical or analytical) to largely reduce the 
errors in Aq^, one could expect an improvement of an order of magnitude in critical 
exponents measures. Nevertheless, we will see below that a truly significant gain will be 
obtained in the L ^ ^ regime. 

A worrying feature of this analysis, is that the scaling corrections for rj and u when 
A = 0.25, seem no longer under control. One would say that we need to simulate in 
significantly larger lattices, to reach the asymptotic regime. We think this to be an effect 
of the competition with the Gaussian fixed-point at A = 0. A similar behavior was found 



in the site diluted Ising model at very small dilution [IS 



6 Numerical Results when ^ <^ L 

As an objective test of scaling, we have chosen the renormalized coupling constant. This 
has been a remarkably difficult quantity to measure until now. Field Theoretic calcula- 



tions [jT4[ yielded 

= 23.73(2), (16) 



which was difficult to reconcile with MC calculations up to correlation-length = 6 



Moreover, calculations |]13| of (see Eqs. (|T^)) show a worrying decreasing tendency, 
even suggesting a hyper-scaling violation. Our data in the right-hand side of Fig. |l| quite 
convincingly show that both the Binder Cumulant and ^/L are non zero and universal 
at criticality. Thus, is certainly not zero, as expected. In fact, using the data and 
procedure of 0, we obtain 

= 5.25(3) (17) 

However, the more interesting quantity to calculate is g^ in the reversed limits ordering. 
The theoretical expectation is that scaling corrections will appear as a power series 

We have simulated with A = 1.0, for correlation-lengths ^ = 1.25,2.5,5 and ^ = 10, 
that correspond to P = 0.316, 0.3588, 0.3729 and 0.3776, respectively. The EMCS consists 
of a Swendsen-Wang update of the spin signs, followed by a Metropohs sweep for changing 
the modulus. Measures were taken every EMCS, using improved estimators (the CPU 
gain using this estimators is not smaller than a factor of 100, if the need for (3 reweighting 
is not as critical as in the FSS region). We have carried out 1.5 x 10® EMCS for each lattice 
size and /3 value. In all cases, but ^ = 10, we have simulated lattices with L/^ ^ 6, 10 and 
12 (we have used lattices L = 8, 12, 16,24,32,48 or 64, depending on ,^). The finite-size 
corrections are expected to drop exponentially with L if periodic boundary conditions are 
used. At our accuracy level (2% in the renormalized coupling constant), no significant 
finite-size corrections are found for L/^ > 6 (see Fig. Thus, for ^ = 10, we have only 
simulated a L = 64 lattice. 



Our data are plotted together with those for the Ising model of Ref. [12 , which have 



been taken at L/^ = 10, thus being quite asymptotic. The horizontal dashed line is the 
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Figure 3: Renormalized coupling constant for (3 < (3c- The Ising data are taken from 
Ref. [0 and correspond to L/C, = 10. The A = 1 data are obtained at L/C, = 6 (dia- 
monds), L/C = 10 (squares) and L/C = 12 (circles). The horizontal dotted line correspond 
to the theoretical value. 



continuum-limit value of Eq. (|l^). For the Ising model we draw a quadratic fit in C~'^ 
(with u = 0.87), constrained to pass through the value in Eq. (^) (horizontal dotted 
line), including only data with ^ > 1. The Ising data at ^ = 0.8, is somehow far from the 
quadratic fit, but for ^ < 1 it is not really surprising that one needs several terms of the 
C''^ series. In the fit, a significantly non zero coefficient for ^""^ is found, as one would 
expect, and also an important quadratic contribution. 

Regarding the A = 1 data, at ^ = 1.25, the scaling-corrections are at the 10% level (to 
be compared with a 50% in the Ising model). From the data we observe a cancellation 
(within errors) of the linear term and a significant reduction of the quadratic corrections. 
This could be interpreted as a full cancellation of the ^""^ series (full decoupling of the first 
irrelevant operator). As for C > 2.5 our results are compatible with the Field-Theoretical 
value, we cannot distinguish if the deviations at ^ = 1.25 are due to other irrelevant 
operators or to the analytical corrections (which would be 0{C~'^^'^))- 

As a final comment, the A = 1.0 two-point correlator, (0(x)0(?/)), does not show a 
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significant improvement on its rotational-invariance when compared with the Ising cor- 
relator. That is not really surprising, because we have included only a first-neighbors 
coupling. So, at correlation-length ,^ ^ 1 we only have a statistical system with highly 
anisotropic couplings, which is naturally reflected in the high-momentum behavior of the 
propagator. 



7 Conclusions 

We have shown that by eliminating the scaling corrections of the scaling-f unction of the 
Binder cumulant in the FSS regime, we obtain a radical improvement of the scaling 
behavior of other critical quantities, both in the FSS regime and in the thermodynamical 
limit. This has been done by considering objective properties of the system under study, 
not depending on an ad hoc RG transformation. This result poses at least three questions. 
How critical is the choice of the tunable parameter. A? (that is, could we have obtained 
such a good result by considering a next-to-nearest neighbors coupling, for instance?) 
How can this optimization strategy be implemented in an asymptotically-free theory? 
(this could be investigated in the two dimensional, 0(3) Non-Linear a model, and will 
be considered in a near future) And last, but not least, will this strategy work in an 
asymptotically-free lattice-gauge theory? 

Regarding the reduction of the errors, the optimal value of A is only approximately 
known, thus the 0{L~'^) corrections cannot be completely disregarded. As a consequence, 
not true gain in the critical exponents determination is achieved. The situation is quite 
better in the thermodynamical regime (L ^ ^), where continuum-results (at the 2% level) 
for the renormalized coupling constant can be obtained at correlation-length ^ = 2.5. Of 



course one could object that a judicious use of the FSS method p6| allows to obtain 
thermodynamical data at very large correlation-length. But also in this case there is a 
real danger that scaling- correct ions will spoil the extrapolation, so a drastic corrections 
reduction will be of benefit. 
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